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Abstract 

The beamforming techniques have been recently studied as possible enablers for underlay spectrum 
sharing. The existing beamforming techniques have several common limitations; they are usually system 
model specific, cannot operate with arbitrary number of transmit/receive antennas, and cannot serve 
arbitrary number of users. Moreover, the beamforming techniques for underlay spectrum sharing do 
not consider the interference originating from the incumbent primary system. This work extends the 
common underlay sharing model by incorporating the interference originating from the incumbent system 
into generic combined beamforming design that can be applied on interference, broadcast or multiple 
access channels. The paper proposes two novel multiuser beamforming algorithms for user fairness 
and sum rate maximization, utilizing newly derived convex optimization problems for transmit and 
receive beamformers calculation in a recursive optimization. Both beamforming algorithms provide 
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efficient operation for the interference, broadcast and multiple access channels, as well as for arbitrary 
number of antennas and secondary users in the system. Furthermore, the paper proposes a successive 
transmit/receive optimization approach that reduces the computational complexity of the proposed 
recursive algorithms. The results show that the proposed complexity reduction significantly improves 
the convergence rates and can facilitate their operation in scenarios which require agile beamformers 
computation. 


Index Terms 

Multiuser MIMO, underlay sharing, generic coordinated beamforming, sum rate maximization, 
fairness, recursive and successive optimization, interference, broadcast, multiple access channels. 

I. Introduction 

The proliferation of wireless devices and ever increasing use of mobile data makes efficient use 
of radio spectrum a critical issue. Cognitive Radio (CR) is identified as a promising technology 
that can facilitate secondary (non-licensed) users to re-utilize or share the licensed bands [HI 
and significantly improves the spectrum utilization. Employing multiple transmit and receive 
antennas at the secondary user (SU) systems can guarantee high spectrum efficiency while 
avoiding interfering with the primary (incumbent) users (PUs). The use of multiple antennas 
offers additional degrees of freedom due to the spatial dimension [|2|-[|4|, which can support 
simultaneous and transparent operation of the SU systems with respect to the PUs (underlay 
sharing). Moreover, multiple antennas provide possibilities for more flexible and efficient resource 
allocation compared to the conventional SU systems (e.g. interweave-based SU systems). 

Recent research activities in spectrum sharing have shown considerable interest in the design 
of practical and efficient beamforming techniques for the underlay-based SU systems [|5l-[[8l. 
The objective of beamforming in the context of underlay spectrum sharing is to maximize the 
SU rate with a given transmit power budget, while keeping the harmful interference to the PU 
system below a predefined threshold. 
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The previous works on beamforming with underlay speetrum sharing generally focus on 
scenarios where only one SU (communication pair or data link) shares the spectrum with one or 
multiple PUs |I51, ll9l- lfT^ . More recently, the design of the beamforming techniques has been 
extended and applied to the multiuser SU scenario [|7]|, |[8l, ifTO . [IT3l - lf2^ . However, there is a 
lack of generic and robust transmit/receive beamforming designs applicable to different system 
models and operating with arbitrary number of antennas and users. The existing beamforming 
techniques are usually designed for a specific system model, i.e. optimized to be used for the 
interference channel IfT^ . [|T9l broadcast channel [f20l . [|^ or the multiple access channel [|2T]| . 
[|22l . Another drawback of the existing schemes is that the maximum number of SUs in the 
system cannot exceed the number of either transmit or receive antennas in the system [|T^ - [|^ . 
Most of literature considers system models and beamforming techniques that are developed for 
the broadcast channel [|71|, IfTOl . lfT3l - lfT5l and the MISO transmit beamforming scenarios. To 
the best of the authors knowledge only Liu & Dong in IfT^ and Scutari & Palomar in IfTSi have 
tried to address the issue of the multiuser interference channel for the MIMO-based underlay 
spectrum sharing. However, the proposed solutions in [|T^ . fITSi can only serve a specific number 
of SUs, which is upper bounded by the number of transmit antennas. 

Even in conventional multiuser MIMO beamforming research, the authors have struggled to 
produce generic multiuser coordinated beamforming techniques relying on convex problems for 
the sum rate optimization. Most of the current work focused on the zero forcing transmit beam¬ 
forming and the optimal maximum ratio combining for the receiver beamforming. Employing 
zero forcing flM . If25l at the transmit receiver side imposes limitations on the number of operating 
users and the number of antennas in the system. Instead of zero forcing, the sum rate optimal 
transmit beamformers for the multiuser MIMO broadcast and medium access channel can be also 
solved using the generalized BC-MAC duality approach [|2^. However, the BC-MAC duality is 
not feasible to underlay spectrum sharing scenarios. The sum rate optimization has been also 
tackled for the MISO interference channel, with the optimal transmit beamformers attained via 
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SINK feasibility and outer polyblock approximation approaches GTl . These approaches result 
in high computational complexity and generally slow convergence. 

In addition, most of the beamforming related underlay sharing literature omits the reverse 
interference, i.e. the interferences caused from the PU to the CR system, in their system models 
and the corresponding beamforming designs Il5l-|[8l, [fT^ - [l2^ . ll^ . In realistic underlay sharing 
scenarios, the reverse interference can severely diminish the performance of the SU system and 
should not be neglected when performing the beamforming optimization 10. In order to deal 
with and align with the PU interference caused to the secondary system, multiple antennas are 
a necessity at the SU receiver sides. 

In this paper we address some of the limitations of the existing work. Specifically, we propose 
several underlay spectrum sharing techniques with respect to the beamforming optimization. 
In particular, the contributions of this paper are the following. First, we propose a generic 
beamforming design by exploiting the advantages of the coordinated beamforming that can be 
utilized for underlay spectrum sharing. Our generic beamforming design can be applied for all 
three major system models, i.e. interference, broadcast and multiple access channels. It allows 
underlay operation of the arbitrary number of users and antennas. Furthermore, the proposed 
beamforming design can be applied to conventional MIMO beamforming systems. Second, 
we propose two novel beamforming algorithms that provide convex solutions for recursive 
computation of the optimal transmit and receive beamformers. The former leverages fairness 
between the users in the system, while the later one maximizes the sum rate in the system. In 
addition, the paper proposes a successive optimization approach that decreases the computational 
complexity of the presented fairness and sum rate maximization beamforming algorithms. These 
algorithms are also applicable to conventional MIMO beamforming systems. Finally, we present 
and employ a more generic system model for the underlay spectrum sharing that takes into 
account the reverse interference caused by the PU system to the SU system. 

The remainder of the paper is organized as follows. Section |II] presents the notation and 
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Fig. 1. System model for multiuser MIMO CR (a) interference channel, (b) broadcast channel and (c) MAC channel 


system model under eonsideration. Seetion |nl] presents the design of the reeursive fairness and 
sum rate maximization algorithms, as well as convexity reduction for successive transmit/receive 
beamforming optimization. Section |IV] assesses the performance of the algorithms under different 
system models considerations and scenario setups. Finally, Section |V] concludes the paper. 


II. System Model 

We consider a multiuser beamforming communication of Ng secondary users operating in 
underlay manner in the same band as Np primary user pairs. The paper considers three possible 
MIMO system models for the secondary user communication, i.e. a multiuser interference (IC) 
channel (Fig. [Tal), a broadcast (BC) channel (Fig. [Tb]) and a multiple access (MAC) channel 
(Fig. [T^. The following text first considers the multiuser MIMO CR interference channel model 
as the most comprehensive one and afterwards presents the respective simplifications for the 
remaining system models. 

The secondary users are equipped with Nt antennas at the transmitter side and antennas 
at the receiver side, while the primary users target a SISO communication. All the channels in 
the system are assumed to follow a CJ\f distribution. In the remaining text, |||| denotes a vector 
norm, tr denotes a matrix trace operation, rank represents a matrix rank operation and blkdiag 
denotes a block diagonal matrix. 
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A. Multiuser MIMO CR interference channel 

In this subsection we consider the most comprehensive MIMO CR interference channel model 
(Fig.fTa]). The secondary system is composed of Ns independent communication pairs (transmitter 
and receiver) aiming to transparently utilize the same band used by Np communication pairs on a 
primary basis. The notations in the subsequent analysis use per communication pair (SU or PU) 
indexing. In particular, the 1x1 channel between the transmitter in the PU pair i and the receiver 
in the PU pair j is denoted as hpp(^ij), where i,j G {l,...,iVp}. The NrXNt communication 
channels between the secondary users are denoted as matrices where k denotes the 

index of the SU pair the transmitter belongs to, I denotes the index of the receiver’s SU pair, and 
k, I E {1,..., Ns}. The NrXl interference channel between the PU transmitter i and SU receiver 
I is represented as hps(i^i), where i E {l,...,iVp} and I E {l,...,iV 5 }, while the interference 
channel between the SU transmitter k and the PU receiver j is represented with the IxA^^ vector 
hsp(fcy), where k E {l,...,iV 5 } and j E {l,...,iVp}. The iV)xl transmit and the N^yil receive 
beamforming vectors for the SU pair k are denoted as rrifc and w^, respectively. 

1) Secondary system perspective: The signal received at the SU receiver I comprises the 
signal from the SU transmitter I summarized with the cumulative interference coming from all 
the PU transmitters, the cumulative interference from the remaining SU transmitters (k 1), 
multiplied by the respective channels, as well as the additive Cff noise at the SU receiver I, i.e. 



k=l;k^l 


where Sp{i),i E 1, ...,Np and Ss{k), k E I, Ns are the symbols transmitted by the PU and SU 
transmitters and ns{i) is the noise at the SU receiver 1. The SU performance can be quantified by 
the Signal-to-Noise-plus-Interference-Ratio (SINR) at receiver side. Considering ([T]), the SINR 
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of the communication of the SU pair I, independent to the symbols realization, is calculated as 


SINRsd) = 


1 


r 

S^Np 1 

Z^i=l 1 


|2 , y^ZV, 1 

1 ^ Z^k=l;k^l 1 

1 w 

1 1 

|2 , 2 

1 +^.(0 


( 2 ) 


where denotes the noise variance at the receiver in the SU pair 1. The SINK at the SU 
receiver I can be also rewritten in the following form 


SINRs^i) = 


Eii + Ek=i;k^i + a; 

Let us define the following NtxNt positive semidefinite matrices 


2 

s(l) 


(3) 


Mfc = mfcinf ^ 0,rank(Mfc) = l,\/k e 
Gssik,i) = y 0, V/c, I e {1, Ns} , (4) 

Gsp{k,j) — ^sp{k,j)^sp{k,j) y 0, rank(Gsp(fcj)) = l,V/c G {1, ...jiV,} ,Vj G {1, ...,Np} . 
Substituting these matrix terms in ([3]), the SINK at the SU receiver I can be calculated as 


SINRs^i) = 


^fc=l;A: 7 ^Z fr(Gss(fc,Z)Mfc) + Ips(l) + 


(5) 


where Ips{i) = E^=i ^fs(i is the summary PU interference received at the SU 

receiver I side. We further define the iV^iV^xiVtiVs block diagonal matrices 


M = blkdiag{Mi,, rank(M) = Ns, 

Xs(i) = blkdiag {Gss(i,z), Gss(i-i,i), [OlATtxATt ’ ^ss(i+i,i), Gss{Ns,i)} J ^ ■■■, ^s} , 

( 6 ) 

Qs(Z) = blkdiag {[OjjVtxAft > •••) [OlxTixiVt > [OjjVjxZVt > •••’ MzVtxAft) > ^ ^ {1; •••) , 

^p(i) blkdiag |Gsp(iy), Gsp(zVgj)]‘ ,j G {1, Np} , 

which conserve the positive semidefinite property of the comprising NtxNt matrices. Here, 
[0]jY^xzvt denotes a zero matrix of size iV^xiV^. 
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Based on the bloek diagonal matriees defined in the SINK at the SU reeeiver I ean be 
simplified to be 


SINRsii) 


tr(Q^(;)M) 


(7) 


All the SU pairs now have a SINK dependenee on the same matrix term M, whieh is a 
rank Ng matrix, eomprising all transmit beamforming veetors information. This derivation will 
prove to have a signifieant impaet onto the eonvexity and resolvability of the optimal transmit 
beamforming veetors. Finally, the power eonstraints of the transmit and reeeived beamformers 
ean be defined as 


m;||' = tr(MO<Pt,V/G{l,...,iV,}, 


( 8 ) 


A typieal eharaeteristie in the multiuser MIMO IC ehannel is the power eonstraint per SU 
transmitter (< Pt), while the reeeiver beamformer is eonstrained to have a max. of unit norm. 

2 ) Primary system perspective: The signal reeeived at the PU reeeiver j eomprises the signal 
eoming from the PU transmitter j, the eumulative interferenee from the all other PU transmitters 
{i 7 ^ j), the eumulative (a.k.a. aggregate) interferenee eaused by all SU transmitters, multiplied 
by the respeetive ehannels, as well as the additive CJ\f noise, i.e. 

Np Ns 

ypU) ~ ^ppUP^pU) + ^ppiip^pU) + +rip(j), (9) 

k=l 

where is the noise at the PU reeeiver j. Sinee the aim of the targeted algorithms in this 
paper is to operate in an underlay manner with the PU system, one of the metries of interest is 
the eumulative SU interferenee power eaused to the PU reeeiver j (irrespeetive to the symbols’ 
realization at the SU transmitter side). Utilizing the matrix terms from (jb]), we define the SU 
interferenee temperature eonstraint for all primary system reeeivers as 

Ns 

^sp{j) ^ ^ 11 11 tr(XpQ'^]\/[) ^ ^ . (10) 

k=l 
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This constraint limits the aggregate interferenee eaused by all SU transmissions (V/ G {1,A^^}) 
to all primary reeeivers (Vj G below a predefined threshold 7 . This is a generie 

underlay speetrum sharing primary system proteetion eonstraint. This subseetion presented 
the respeetive derivations (SUs SINK and PU interferenee eonstraint) for the multiuser MIMO 
CR interferenee ehannel. Considering Fig. [Tal it is intuitive that there is no limitation on the 
number of operating seeondary links Ng, i.e. supported eoncurrent data streams. The eombined 
virtual MIMO system has NgXNt transmit and NgXNr reeeive antennas, meaning that the degrees 
of freedom in this system ranges with the number of SU pairs. 

B. Multiuser MIMO CR broadcast and multiple-access channels 

The multiuser MIMO broadeast and multiple-aeeess ehannels are analogous to a seenario of 
single eell MIMO transmissions. A broadeast ehannel is the ease of single base station (BS) 
serving multiple terminals in downlink (Fig. (Ta]), while a MAC channel is the case of multiple 
terminals eommunieating with a single BS in uplink (Fig.fTbl). The following subseetions present 
the required simplifieations to adapt the approaeh to these system models. 

1) BC channel: In the ease of the broadeast ehannel (Fig. [Tbl) . the SU system is eomposed of a 
signal souree (transmitter) and multiple signal destinations (reeeivers). In this ease, the SU signal 
and the SU interferenee to a single SU receiver arrive via the same ehannel, multiplied with 
the respeetive transmit beamforming veetors. The similar logie applies to the SU interferenee 
eaused to the PU reeeivers. The following equations represent the differenees (simplifieations) 
with respeet to the multiuser MIMO interferenee ehannel (eqs. dH) and ®) 

Gss(i,/) = Gs5(2,i) = ... = Gss(7v^^Z), V/ G {!,... 

( 11 ) 

^sp(2,i) ••• hsp(Afs7)5 Vj G {1,..., , 

Gsp(ij) G^p( 2 j) ••• G^p^jVsj)5 V) ^ 
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For the multiuser MIMO broadcast channel, the transmit power constraint in ([8]), is defined as 

Ns 

5^||mz||'=tr(M)<Pi (12) 

1=1 

As can be noticed, the power is shared among the transmit beamformers (i.e. SU data streams), 
providing options for transmit power allocation optimization. Finally, the broadcast channel 
imposes a limitation of the number of served secondary data streams (Ng). Considering equation 
(fTTI) and Fig. [TQ it is intuitive that the degrees of freedom are limited by the number of transmit 
antennas Nt. The combined virtual MIMO system has Nt transmit and NgxNr receive antennas, 
meaning that the number of served SU data streams must be Ng < Nt. 

2) MAC channel: In the case of the MAC channel (Fig. [1^, the SU system is composed 
of multiple signal sources (transmitters) and single signal destination (receiver). All SU signals 
and SU and PU interference arrive at the same receiver, imposing the following adaptations and 
simplifications with respect to the multiuser MIMO interference channel (eqs. dH) and db])) 


Hss(fc,i) Hss(fc,2) ••• , Vfc G ■( 1,..., } , 

••• ^ps(i,Ns) ^ •■■yNp^ . 


(13) 


The transmit power constraint is the same as in case of the MIMO interference channel. Again, 
the MAC channel imposes a limitation of the number of served secondary data streams (Ng). 
Considering Equation (fTSI) and Fig. [T^ it is clear that the degrees of freedom are limited by the 
number of receive antennas N^.. The combined virtual MIMO system has NgXNt transmit and 
Nr receive antennas, implying that the degrees of freedom in this system is limited by Nr. This 
means that the number of independent and concurrent secondary data streams must be Ng < Nr. 


C. Summary 

In general, the same SINR related equations d2l)-d2]) and PU protection related expressions 
(fTOl) can be used for all three models, with the respective simplifications shown in (fTTI) and (fT3]) 
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for the multiuser BC and MAC models, respeetively. Finally, a generic transmit power constraint 
for all considered channel models is defined as 


Ns 


limin'= tr(M0 <Pt,V/e ,iV,} or V ||mzf = tr(M) < 


(14) 


IC and MAC channels 


1=1 


BC channel 


and the generic receive power constraint as 




(15) 


Besides the previously denoted limitations on the number of served secondary data streams, all 
system models have a limitation on the number of protected PU receivers. In particular, in order 
to provide a secondary communication, the maximum number of protected primary receivers is 
limited by Nt — 1, for each of the three multiuser MIMO system models, irrespective on the 
number of SU transmitters. This limitation is caused by the interference alignment capabilities, 
i.e. such system is able to provide transmit beamforming based interference alignment of up to 
maximum of Nt streams. Note that this limitation, along with the number of served SU data 
streams limitations are all system model limitations. All simplifications, limitations and scenario 
specifics are summarized in Table HI Further adaptations can be easily made to cover additional 
scenarios (system models), such as the multicell multiuser MIMO scenario. In the next section 
we will propose two generic multiuser beamforming strategies for underlay spectrum sharing. 


TABLE I 

Multiuser MIMO system model specifics for coordinated beamforming optimization 


System model 

Simphfications 

Limitations 

Capabilities 

Interference channel 

/ 

Np<Nt-l 

transmit/receive bf. 

Broadcast channel 

in (1111 

Np<Nt-l 

Ns<Nt 

transmit/receive bf. 
power allocation 

MAC channel 

in (1131) 

Np<Nt-l 

Ns < Nr 

transmit/receive bf. 
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III. Generic multiuser beamforming algorithms for underlay sharing 

In this section we utilize our results from the previous section in order to design two generic 
multiuser beamforming algorithms for underlay spectrum sharing. In particular, the algorithms 
assume operation in a centralized manner, i.e. cooperation between the SU nodes. All channels 
are assumed to be frequency flat, fixed, and perfectly known (full channel state information) 
prior to the optimization. This is a justified assumption, since each of the channels can be 
estimated prior to the optimization and the operation. In this section we specifically will derive 
a fairness optimization algorithm where the radio resources are shared among SUs to achieve 
balanced SINK distribution per secondary data stream and a sum rate maximization algorithm 
that maximizes the sum rate of the secondary system. 

Since the problem of finding the optimal transmit and receive beamforming vectors simulta¬ 
neously is a non-convex problem, both proposed algorithms treat the optimization iteratively. In 
particular, the general idea is to divide the transmit and the receive beamforming optimization 
into two separate sub-problems and optimize the beamforming vectors recursively, i.e. calculate 
the transmit beamformers for known (previously optimized) receive beamformers and vice versa. 


A. Multiuser CR beamforming fairness algorithm 

The multiuser beamforming fairness optimization for underlay spectrum sharing can be real¬ 
ized using SINK max-min optimization with respect to the primary system and power constraints. 
Utilizing the SINK expression in (|3]), the optimization problem can be formulated as 


{m/, = arg max min [SINRs(i)} , 

subject to: (fTOl) . (fT4l) and (fT5l) . 


(16) 


In particular, this optimization problem tends to balance the SINK distribution per SU data 
stream, keeping the aggregate interference caused to all the primary receivers below a predefined 
threshold 7 , and the transmit power of SU transmitter(s) below a maximum of Pt. The receive 
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beamformers have a unit norm constraint. However, the simultaneous optimization of both the 
receive and transmit beamformers, results in a non-convex optimization problem. In order to 
transform the problem into a convex one, we divide the optimization problem into two sub¬ 
problems, i.e. the transmit beamforming and the receive optimization sub-problems. 


1) Transmit beamforming optimization sub-problem: Assuming that the receive beamforming 
vectors are known, a transmit beamforming fairness optimization sub-problem can be formulated 
using the SINK definition in (|7]) and the respective constraints in (fTOl) and (fT4l) . This yields 


M = arg max min 

M I 


tr(Q^(i)M) 

tr(X,(;)M) + ^ps{l) + 




subject to: rank(Mi) = 1,V/ G {1, ...,Ns} , (fTOl) and (fT4l) . 


(17) 


The problem in (fTTI) can be transformed into a convex one with a semidefinite relaxation, i.e. 
omitting the rank constraints. It can be solved with Dinkelbach fractional programming type 
of algorithms Il29l and using common semidefinite optimization tools. The aim is to find the 
positive semidefinite block diagonal non-full rank NtNgxNtNg matrix M, and extract the optimal 
transmit beamforming vectors accordingly in each iteration of the algorithm. This iteration based 
Dinkelbach algorithm in the n-th iteration operates on the following optimization problem (which 
is equivalent to a semidefinite relaxed problem in (fTVl) ): 


= argmax(r) 

^ (18) 
subject to: tr(Q^(i)M) - 5^i„(tr(X^(i)M) -f > r, ^ and (El), 

where 6min = min ..., | is the minimal SINK acquired in the previous 

77,-1 iteration of the algorithm, and r represents the fairness coefficient, i.e. the difference 
between the maximal and the minimal SINK. The transmit beamformers are calculated 

as the principal singular vectors of the positive semidefinite matrices extracted from 

the block diagonal matrix M obtained by the transmit beamforming optimization in (fTSl) . 
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2 ) Receive beamforming optimization sub-problem: The optimal unit norm receive beamform¬ 
ing vectors that maximize the SINK per SU data stream I (eq. (|7])), for known transmit 

beamforming vectors can be calculated as a Rayleigh quotient: 


W; = 






^ -U /t2 


-1 




(19) 


3) Algorithm for multiuser CR beamforming fairness optimization: The transmit and receive 
optimizations presented in previous subsections, can be programmed as a recursive fairness 
optimization algorithm in the following manner. The optimal beamforming vectors are calculated 
in a cyclic manner, i.e. calculating the transmit beamformers for known receive beamformers, 
and vice versa, calculating receive beamformers for known transmit beamformers until iteratively 
reaching a globally optimal solution. Algorithm [T] presents the pseudo code of this multiuser CR 
beamforming fairness optimization algorithm. 


B. Multiuser CR beamforming sum rate maximization algorithm 


Employing Shannon capacity, the achievable capacities for the SU data stream I and the entire 
secondary system are 


Cs(i) = log2(l + SINRs(i)), 

Ns 

1=1 


( 20 ) 


Based on the SINR expression in ([3]) the beamforming optimization problem for the underlay 
sharing can be defined to maximize the sum rate of the secondary data streams complying with 
the interference constraints, i.e. 


Ns 


{mz, = arg 


max 


+ SINR,(i^) 


1=1 


subject to: (fTOl) . (fT4l) and (fT5l) . 


( 21 ) 
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Algorithm 1 Multiuser CR beamforming fairness optimization 

STEP 1: Initialization 

n ^ I 

for I = 1 : Ns do 

<— random feasible vectors 
■«— based on 0 

STEP 2: Minimum SINR calculation 

n •<— n + 1 

^ min{5/7Vi?J-'\ 5/iVi?} 

STEP 3: Transmit beamforming optimization 

^ based on (fTSl l 
for I = 1 : Ng do 

m|"^ - 1 — principal singular vector of 

STEP 4: Receive beamforming optimization 

for I = 1 : Ns do 

^ based on ( fT9] l 
STEP 5: Stopping criteria 
if r > e then 

for I = 1 : Ns do 

<r- based on 0 
go back to STEP 2 
else 

for I = I : Ns do 

(n) (n) 

^ w/ w; 

Stop algorithm 


In this optimization we aim to maximize the summary seeondary system sum rate, while keeping 
the aggregate interferenee caused to all the primary receivers below a predefined threshold 7 , the 
transmit power of all SU transmitter(s) below a maximum of P* and the receive beamformers 
power below one. The Equation (|2TI) is a non-convex optimization problem. Therefore, we 
consider several techniques to alleviate the non-convexity of the optimization problem and 
transform the problem into a convex one. Similarly to the fairness optimization, we divide 
the problem into two separate sub-problems and optimize the beamforming vectors recursively. 


1) Transmit beamforming optimization sub-problem: In order to transform the problem in 
(I 2 TI) into a convex one, as mentioned, first, the problem is considered only from the transmit 
beamforming optimization perspective, utilizing the SINR expression in (|7]) and performing a 


April 27, 2015 


DRAFT 





16 


semidefinite relaxation 


M 


arg max 

M 


Ns 

5^l0g2 


1=1 



tr(Q^(;)M) 

tr(X,(oM) + Ips{l) + cr 



subject to: (fTOl) and (fT4l) . 


( 22 ) 


The sum rate of the secondary data streams can be represented as difference between two sums 
of logarithmic terms 


C', = /(M)-i7(M), 

Ns 

f (M) = ^ log 2 (tr ((Qs(i) + Xs(z))M) + Ips(i) + o-s(z)) , ^23) 

1=1 

Ns 

t/(M) = ^ logs (tr (Xs(z)M) + /p,(z) + aj;)) . 
z=i 

Both terms, /(M) and are concave functions, and therefore, the optimization problem in 

(l22l) is still a non-convex one. In order to alleviate the non-convexity we use the same approach 
as in If^ . In particular, due to the properties of the logarithm function, (^(M) is weakly sensitive 
to changes of the variable M, so in a local (and fairly large) neighborhood of a random point 
M(”) we can use its linear (first order Taylor) approximation 


c/(M) ^ c/(M(”)) + (Vc/(M(”)), M - , (24) 

where V 5 '(M) is the first order derivative of the function ^'(M) in point calculated as 

V9(M) = ^ ---. (25) 

,=i (tr (X.,,,M) + 4.,,, + 111 2 

Considering (|2^ . (I24l) and (1251) the sum rate function in the local neighborhood of can be 
approximated as 

C, ^ /(M) - c/(M(”)) - (Vc/(M(”)), M - , (26) 
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which is now a concave function and has a single optimum In order to derive a sum 

rate optimization problem that will globally eonverge, we have to eonsider an important property 
of the funetion (y'(M). Sinee this funetion is eoneave, its gradient is also its super gradient, i.e. 
5 f(M) has the following property 

^(M) < M - . (27) 


This property results in the following eorollary 




(28) 


This means that the new solution is always better than the previous Sinee that 

eonstraint set in (l22l) is eompaet, starting from a random feasible point, the optimization will 
always find a better feasible point. The sum rate optimization problem in (l22l) ean be now solved 
iteratively, with the optimal in the n-th iteration ealeulated as an output of the following 
eonvex optimization: 


= argnmx{/(M) - c/(M(”-i)) - (Vc/(M(”-i)), M - , 

subjeet to: (fT^ and (fT^ . 


(29) 


Similar to the fairness algorithm, the transmit beamformers are ealeulated as the 

prineipal singular veetors of the positive semidefinite matriees extraeted from the bloek 

diagonal matrix M obtained by the transmit beamforming optimization in (l22l) . 


2 ) Receive beamforming optimization sub-problem: The unit norm reeeive beamforming vee¬ 
tors that maximize the SINK per SU data stream is ealeulated in the same manner as 

in the fairness algorithm, using formula (fT9l) . 
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3) Algorithm for multiuser CR beamforming fairness optimization: The recursive (trans¬ 
mit/receive beamforming) sum rate optimization algorithm for underlay spectrum sharing is 
defined in the following manner. The optimal beamforming vectors are calculated in a iterative 
manner, i.e. calculating the transmit beamformers for known receive beamformers, and vice 
versa, calculating receive beamformers for known transmit beamformers, until reaching a globally 
optimal solution. Algorithm [2] presents the pseudo code of this multiuser CR beamforming sum 
rate optimization algorithm. 

In order to reduce the complexity of the recursive transmit beamforming optimization and 


Algorithm 2 Multiuser CR beamforming sum rate optimization 

STEP 1: Initialization 

n <— 1 

for I = 1 : N. do 

, Wj" <— random feasible vectors 
SINR^f'^^ <r- based on (|3]l 
•<— based on (|4]i and (|6]l 
■<— based on (l20t 
STEP 2: g, Vg calculation 
n •«— n + 1 

•<— based on (l2?t 
Vg(M(”-')) ^ based on dH 

STEP 3: Transmit beamforming optimization 

^ based on (l29l l 
for I = 1 : Ns do 

m|"^ <r- principal singular vector of 

STEP 4: Receive beamforming optimization 

for I = 1 : Ns do 

^ based on (fT9] l 
STEP 5: Stopping criteria 

for I = 1 : Ns do 

SINR^fj.^ <r- based on @ 

<r- based on (l20l i 

if > e then 

go back to STEP 2 
else 

for I = I : Ns do 

(n) (n) 

; w/ ^ 

Stop algorithm 


the overall complexity of the multiuser CR beamforming underlay sharing (Algorithms [U and 
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[2l), this paper considers an additional simplification of the transmit beamforming sub-problem. 
Applying the generalized Kantorovich and Cauchy Schwartz inequalities, the SINK in eq.©, 
can be easily shown to be lower bounded only by the transmit beamformers, i.e. 


SINRsd) = 


+ Ek=i-,k^i + a 


2 

s{l) 


- v-Af. 




(30) 


™ ,v -= SINRi%. 

Eiii ^?s{i,i)Ksii,i) + Efc=i;Mz 


Considering that S'/A^i?s(z) is lower bounded by SINR^^^^jy we can do the transmit beamforming 
optimization just based on the based on the expressions. Therefore, we define the 

following simplified terms 






^ O.Vfc,; 6 {1,.... iV,} , 

K() = ““‘“g {G?.(U)- . G“Ci,„..... ^ 0. VI e {1, JV.} . 


2 = 1 

fB 


Qs(f) “ blkdiag | , •••, [0] 


NtxNt ’ '^ss(Z,/)’ L^j NtxNt 


JNtxNt 


y0,Vle{l,...,N4. 

(31) 


Now the transmit beamforming optimization can be performed irrespective to the receive beam¬ 
forming, using these terms in the optimization problems provided by (fTSl) and (1281) . The receive 
beamforming optimization is solved successively after the transmit beamforming optimization is 
finished. This will reduce the the number of matrix operations per iteration (refer to Algorithms 
[T] and ©, but will result in sub-optimal behavior. 


C. Summary 

This section presented two multiuser beamforming algorithms for underlay operation in pri¬ 
mary system frequency bands, i.e. fairness and sum rate maximization algorithms, referred to 
as Fairness and SRM in the remaining text, respectively. Both algorithms can be applied to all 
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three MIMO system models, irrespeetive on the number of antennas in the system. Although 
the BC and MAC system models impose limitation on the number of served seeondary data 
streams Ng, the sum rate maximization algorithm can also work with higher number of secondary 
users, optimally selecting the subset of data streams with lowest channel correlation. Thus, this 
algorithm opportunistically maximizes the overall sum rate of the secondary system. In these 
cases, the fairness optimization algorithm will significantly decrease the overall sum rate (the rate 
per user). The following section will present the performance evaluation and prove the benefits 
of the employment of these algorithms in representative scenarios. 


IV. Performance analysis 

In this section we assess the performance of the proposed beamforming algorithms by fo¬ 
cussing on fairness and sum rate maximization. The performance analysis is performed in terms 

Ns Jl 

of the sum rate {Cg = X] log 2 (l + SINRg^i))) and the normalized sum rate (0^ = Cg/Chounda 

i=i 

We investigate also the benefits of exploiting the complexity reduction approach for successive 
optimization for both beamforming algorithms in terms of the Convergence rate. Moreover, 
the section introduces the SNR deviation (i.e. signal deviation, log-normally modeled) as an 
evaluation metric. The SNR deviation is a valuable metric that reflects the realistic behavior of 
the system. In practice it is almost impossible for all or group of secondary data streams to 
experience equal SNRs on their communication channels. The assumption of constant SNR at 
the users has been frequently utilized in previous works and can lead to significantly incorrect 
conclusions. Our scenario parameters are given in Table HIl 


^Cbound denotes the sum rate point-to-point outer bound of the Ns active data streams, which assumes no SU and no PU 
interference in the communication and can he computed as defined in dn. In order to provide easier tractability, for the 
remainder of the paper the term Cbound will be denoted as the upper bound. 
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Parameters 

Values 

Simulation environment 

Matlab 

Optimization tool 

SDP via CVX 

Stopping threshold (e) 

Tir^ 

SNR [dB] 

0 ; 30 

SNR deviation (cr) [dB] 

0 : 5 : 10 

Channel model 

Circular complex normal 

Transmit antennas Nt 

2:2:6 

Receive antennas W 

2:2:6 

PUl threshold (7) [mW] 

Hr™ 

Number of PUs 

1 : 1 : (M - 1) 

Monte Carlo trials per configuration 

1000 


A. Interference channel 


From system model and design perspeetive, the interferenee ehannel represents the most 
generic and complex scenario where all active pairs interfere between each other. In order to 
provide efficient communication over this type of channel, there must exist cooperation between 
the active nodes in the network. From the practical perspective, this scenario can be related to 
coordinated multi-cell communication as well as ad-hoc networks where multiple independent 
and physically collocated transmissions occur at the same time frame. 

Fig. [2^ depicts the achieved sum rate of the secondary system for the fairness and sum 
rate maximization (SRM) beamforming algorithms versus the SNR. We consider also different 
antenna configurations. It is evident that the SRM algorithm always outperforms the Fairness 
algorithm for any SNR and antenna configuration. However, the performance gain of the SRM 
over the Fairness algorithm diminishes for higher number of antennas. This behavior is a result of 
the higher number of available degrees of freedom in the system, facilitating improved sum rate 
performance per user and thus more efficient operation of the Fairness algorithm. In addition. 
Fig. |2^ compares the performance of both algorithms with respect to the upper bound. By 
increasing the number of antennas in the system, i.e. the degrees of freedom, both algorithms 
improve their performances and operate with rates that approach the upper bound. 
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Fig. 2. a) Sum rate vs SNR, for different antenna configurations (Nr = 2, Ns = 2, Np = 1, a = 5dB) b) Sum rate vs 
Number of SU devices, for different SNR values (Nt = 4 Nr = 2, Np = 2, a — OdB) 

Fig. [2b] depicts the sum rate of the secondary system versus the number of active secondary 
users and different SNR values. The SRM algorithm always achieves higher sum rates compared 
to the Fairness algorithm. Moreover, for higher number of secondary users the performance gap 
between the beamforming algorithms increases. This is a result of the specific optimization goal 
of the underlying beamforming process. For the SRM case, the algorithm favors and primarily 
maximizes the sum rate of the data streams with favorable channels conditions (those that 
experience a lower PU and cross SU interference, higher SNR, etc.) neglecting the remaining 
data streams’ rates. In contrast, the Fairness algorithm strives to leverage the performance of 
all users and induce fairness among them, regardless of the attained channels conditions. For 
high number of secondary users the amount of pairs with unfavorable channels conditions will 
increase, thus forcing the Fairness algorithm to improve the sum rates of those pairs at the price 
of decreasing the overall sum rate. 

It is also evident from Fig. [2b| that both algorithms can serve arbitrarily large number of users, 
and are not bounded by the number of transmit and receive antennas. By carefully analyzing 
([7]) and Fig. [Tal it is apparent that the solution space dimension scales with Ng thus, providing 
the possibility to serve any number of secondary users. Nevertheless, the sum rate of both 
algorithms decreases in comparison to the upper bound as Ng increases. This conclusion is 
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counter intuitive with respect to the model presented in Equation ([7]), where the optimization 
process is centralized and exploits the information of all SU channels in the system. However, 
in terms of the transmission strategy for the secondary user streams, the optimization process is 
restricted only to take into consideration and utilize the antennas of the respective transmitter 
receiver pair. This limits the possible solution set for the optimization process and results in 
lower performance when compared with the upper bound. 

B. Broadcast channel 

This subsection evaluates the performance of both algorithms for the BC channel. The BC 
channel represents a specific multiuser scenario realization comprising one transmitter and mul¬ 
tiple receivers. From the practical perspective the most frequent appearance of the BC channel 
is the downlink cell communication. 

Fig. [3^ shows the achieved sum rate of both beamforming algorithms versus the SNR for 
different SNR deviations. As discussed at the beginning of this section, the SNR deviation is a 
valuable metric that can enable more comprehensive analysis for the behavior of the beamforming 
algorithms. Fig. [3^ indicates that the SNR deviation does not impair the achieved sum rate of the 
SRM algorithm. However, it has significant impact on the performance of the Fairness algorithm, 
i.e. higher SNR deviations lead to decreased sum rates. Since the Fairness algorithm strives to 
balance the performance between the data streams, it will exploit the available resources on 
improving the rates of the data streams that experience lower SNR values and ultimately achieve 
lower sum rates compared to the SRM algorithm. 

Fig. [3b] depicts the achieved sum rate of both beamforming algorithms versus the SNR for 
different antenna configurations. Similar conclusions stand as for Fig. [2^ The SRM algorithm 
always outperforms the Fairness algorithm for any SNR and antenna configuration. With the 
increasing the number of antennas in the system, both algorithms improve their performances and 
operate with rates that approach the upper bound. Fig. also shows that the SRM performance 
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Fig. 3. (a) Sum rate vs SNR for different SNR deviations (Nt = 6 Nr = 4, Ns = 2, Np = 2), (b) Sum rate vs SNR for 
different antenna configurations (Nr = Na = 2, Np = 1, a = lOdB), (c) Normalized sum rate vs Number of PUs, for 
different SNR deviations (Nt = 6 A^r = 4, = 4, SNR = lOdB) 


gain (the sum rate differenee) over the Fairness algorithm, is independent of the number of 
antennas. This is opposite to the eonelusions dedueted from the IC ehannel seenario. For the 
BC ehannel there exists a single/eentralized transmission point that provides the possibility to 
alloeate the available degrees of freedom and transmit power per user in a more effieient manner. 
This enables the SRM algorithm to opportunistically adapt the transmit beamformer gain for 
eaeh user in terms of the obtained ehannel eonditions, and exploit the transmit power dimension 
in order to maximize the sum rate. 

Fig. [3c] depicts the normalized sum rate 0^ of both algorithms in relation to the number of 
active PUs and different SNR deviations. As seen from the figure, both algorithms achieve lower 
sum rates for higher number of PUs. The increased number of PUs, stipulates the beamforming to 
exploit most of the available degrees of freedom to mitigate the inter PU-SU system interference 
and lead to decreased SU sum rates. However, for a small number of PUs both algorithms 
operate with rates that are very near the upper bound. With regards of the SNR deviation both 
algorithms exhibit contrasting behaviors. This is due to the reverse interference. For a small 
number of PUs the scenario is dominated primarily by the SNR and the same conclusions hold 
as for Fig. [3^ In contrast, for a large number of PUs the scenario is dominated by the reverse 
interference caused by the PU systems. Since the SRM algorithm predominantly maximizes the 
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Fig. 4. (a) Sum rate vs SNR for different SNR deviations (Nt = 6 Nr = 4:, Ns = 2, Np = 2), (b) Sum rate vs SNR for 

different antenna configurations (Nr = 4, Na = 2, Np = 1, a = lOdB), (c) Normalized sum rate vs Number of PUs, for 
different SNR deviations (Nt = 6 Nr = 4, Ng = 4, SNR = lOdB) 

rates of the data streams with better ehannel eonditions (in this ease higher SINK), bigger SNR 
deviations will result in improved opportunities for the rate maximization process. With respect 
to the Fairness algorithm, the interference limited operation alleviates the SNR deviation effect 
on the algorithm’s performance, as the majority of the SUs will exhibit unfavorable channels 
conditions (low SINRs) and result in convergence of the sum rates. Fig. |3c] also shows that both 
algorithms can be applied to a conventional BC channel beamforming scenario, i.e. BC channel 
without PU systems (Np = 0). 

C. Multiple access channel 

In this subsection we consider the performance of both algorithms for the MAC channel. The 
MAC channel represents a specific multiuser scenario realization comprising multiple transmitters 
and one common receiver. From the practical perspective, the most frequent appearance of the 
MAC channel is the uplink of cell communication systems. 

Fig. 11^ depicts the achieved sum rate of both beamforming algorithms versus the SNR and 
different SNR deviations. The same conclusions stand as for Fig. [3^ More specifically, the SNR 
deviation does not impair the achieved sum rate of the SRM algorithm, however, for the Fairness 
algorithm, higher SNR deviations lead to decreased sum rates. 

Fig. |4b] depicts the achieved sum rate of both beamforming algorithms versus the SNR and 
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different antenna eonfigurations. The SRM algorithm always outperforms the Fairness algorithm 
for any SNR and antenna configuration. With the increase of the number of antennas in the 
system, both algorithms improve their performances and operate with rates that approach the 
upper bound. Additionally, Fig. |4b] shows that the SRM’s performance gain over the Fairness 
algorithm diminishes for a higher number of antennas as a consequence of the distributed 
transmission system model (see, discussions for Figj^ and Fig. |3bl) 

Fig. |4^ depicts the normalized sum rate 0^ of both algorithms as a function of the number of 
active PUs and different SNR deviations. The same conclusions stand as for Fig. [3^ Specifically, 
when the number of PUs is increased both algorithms achieve lower sum rates. Moreover, due 
to the interference limited operation of the SU system, higher SNR deviations improve the SRM 
algorithm’s sum rate, but however, have almost no effect on the Fairness algorithm. For a small 
number of PUs both algorithms operate with rates that approach the upper bound. As a result of 
the SNR dominated operation of the SU system, higher SNR deviations decrease the performance 
of the Fairness algorithm, while having negligible impact on the SRM’s performance. Fig. |4^ 
also shows that both algorithms can be applied to a conventional MAC channel beamforming 
scenario, i.e. MAC channel without PU systems (Np = 0). 

The presented results in Sec.IV.A, Sec.IV.B and Sec.IV.C, show that both algorithms are 
capable of efficient operation in all three system models (i.e. IC, BC and MAC), hence verifying 
the generic applicability of the proposed beamforming design (Sec. HU). Moreover, by comparing 
the results from Sec.IV.B and Sec.IV.C, it can be concluded that the SRM and the Fairness 
algorithms operate equally efficient in both downlink (BC channel) and uplink (MAC channel), 
demonstrating their universality. With respect to the achieved sum rates, the SRM algorithm 
always outperforms the Fairness algorithm, regardless of the system model and scenario setup. 
However, this is an anticipated result as the Fairness algorithm focuses on leveraging the rates 
between all data streams at the price of decreasing the overall sum rate. 
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Fig. 5. Convergence rate comparison (Nt = 4 Nr = 4, Ns = 10, Np = 2, SNR = lOdB, a = 5dB) 
D. Complexity reduction 


This subsection assesses the convergence behavior of the SRM and the Fairness optimization, 
for the recursive and successive optimization cases. It verifies the advantages of the complexity 
reduction approach (Sec.III.C) in terms of the Convergence rate of both algorithms. For the 
remainder of this subsection the successive SRM and Fairness optimization algorithms will 
be denoted as SRMRed and FaimessRed, respectively. Fig. |5] depicts the Convergence rate of 
the SRM, Fairness, SRMRed and FaimessRed algorithms, for the IC channel scenario. The 
convergence rate is presented as the residual error (the difference between the attained value 
and the optimal valued) in function of the number of iterations performed by the beamforming 
algorithm. The IC channel is specifically chosen for the performance evaluation in Fig. [51 as the 
most complex of the three system models, and most suitable for emphasizing the convergence 
behavior of the beamforming algorithms. 

Fig. |5] shows that the SRM algorithm requires a higher number of iterations in order to 
converge to the optimal solution when compared to the Fairness algorithm. It is also evident 
that the complexity reduction approach significantly improves the convergence rates of both 
beamforming algorithms, and can facilitate their operation in scenarios which require agile 


'The values of interest represent the optimization functions values for the SMR l l29t and the Fairness l llSt optimization, i.e. 
the sumrate and the fairness coefficient r, respectively. 
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beamformer computation. Additionally, Fig. [5] shows that the SRMRed algorithm achieves higher 
eonvergenee rate in eomparison to the FaimessRed algorithm. 

This section evaluated the performanees of the SRM and the Fairness algorithms. It verified 
their universal applieability and rate effieieney for for variety of parameters and seenario se¬ 
tups. Moreover, it validated their underlying eomplexity and possible reduetion by applying the 
approaeh presented in See.III.C. 


V. Conclusion 

Reeently, the eoneept of MIMO and beamforming has been extensively investigated as a 
possible enabler of the underlay spectrum sharing. A common limitation of earlier proposed 
beamforming teehniques is that they either eannot exploit arbitrary number of antennas in 
the system and/or eannot serve an arbitrary number of seeondary data streams in the system. 
Moreover, the existing beamforming teehniques frequently treat the reverse interferenee as noise 
at the seeondary reeeiver. Additionally, all of the existing beamforming works we are aware 
of that focus on underlay speetrum sharing propose algorithms that laek universality and are 
speeifieally designed either for the interferenee ehannel, broadeast ehannel or the multiple aeeess 
channel. This segregated design signifieantly limits their applieability in praetieal seenarios, 
where the underlying beamforming proeess should be eapable of operating in any possible 
multiuser seenario. 

In this paper we have proposed several advancements and novel techniques with respect to the 
beamforming based underlay speetrum sharing. First, the paper extends the eommonly utilized 
underlay spectrum sharing model by ineorporating the reverse interferenee in the beamforming 
proeess. Additionally it proposes a generie eombined beamforming design that is applieable for 
any multiuser seenario and is applicable for eonventional as well as for underlay speetrum sharing 
based systems. Furthermore, the paper develops two novel multiuser beamforming reeursive 
algorithms for user fairness and sum rate maximization based on the proposed generie beam- 
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forming design. Both algorithms provide operate on convex problems for the computation of the 
optimal transmit and receive beamformers. The paper also elaborates on a possible successive 
optimization approach that decreases the complexity of the proposed fairness and sum rate 
maximization beamforming algorithms. The presented results in the paper clearly show that 
both beamforming algorithms are capable of efficient operation in any multiuser scenario. Both 
algorithms can operate for any number of antennas and users in the system, thus verifying their 
universality and practical applicability. Moreover, the results show that the proposed complexity 
reduction approach significantly improves the convergence rates of both algorithms and can 
facilitate their operation in scenarios which require agile beamformer computation. 
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